Hadron yield correlation and constituent quark degree of freedom in heavy ion collisions 
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Based on the assumption of the production of deconfined quark matter, we use a quark combination model to 
systematically investigate hadron yields in heavy ion collisions from RHIC ^fs^N - 200, 130, 62.4 GeV to SPS 
Ebeam = 158, 80, 40, 30, 20 AGeV. We find that as the collision energy is greater than or equal to 30 AGeV the 
yields of various hadrons, their correlations, in particular, the observables A = jj-rj, and B = ^ t_ , are all 
reproduced; however, as the collision energy drops to 20 AGeV quark combination fails. This indicafes that the 
constituent quark degrees of freedom represent a decisive factor in thermal hadron production above 30 AGeV 
and seem to be invalid at 20 AGeV. In addition, hadron yields as well as particle ratios at midrapidity in the most 
central Pb+Pb collisions at \/snn =5.5 TeV are predicted. 
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I. INTRODUCTION 



The question at which collision energy in heavy ion 
collisions the deconfinement is first reached has attracted 
more and more attentions in recent years fl-Sl- The 
Beam Energy Scan programme of NA49 experiment at 
the CERN-SPS has suggested a preliminary answer — 
around 30 AGeV|H. The ongoing Beam Energy Scan 
programme of STAR Collaboration at Brookhaven Na- 
tional Lab provides an opportunity to study it in more 
detail lfiol [Till . Once the deconfined hot and dense quark 
matter is produced in heavy ion collisions, the observables 
of various thermal hadrons after hadronization, e.g. yields 
and momentum spectra etc, have some coiTelations orig- 
inated from early quark degrees of freedom. One of the 
most typical examples is the elliptic flow (v2) of hadrons 
measured at RHIC energies. As both v2 and transverse 
momentum(p7) are divided by the constituent quark num- 
ber of hadron, the rescaled v2 of various baryons and 
mesons, which are just that of constituent quarks, almost 
coincide with each other in the intermediate pj range lfl2ll . 
If the hot and dense quark matter is hadronized by quark 
(re-)combination/coalescence iT3l - [l5ll . as is commonly ac- 
cepted, these correlations of hadrons can be beautifully ex- 
plained. In quark (re-)combination/coalescence scenario, 
quarks and antiquarks are available in unbound state before 
hadronization and they can coalesce freely into various 
hadrons, and thereby these correlations from early quark 
degrees of freedom among different hadron species are nat- 
urally formed. On the other hand, if the deconfined quark 
matter is not produced at all in collisions, there is no free 
quarks and antiquarks (much less their subsequent com- 
bination) and these so-called "quark-level" correlations of 
hadrons maybe disappear or contort. Therefore, we can 
study whether the deconfinement is achieved by investi- 
gating these correlations among various hadrons produced 
in heavy ion collisions. 



Hadron yield is one of the most significant observables 
from which one can obtain a lot of important information 
on the hot nuclear matter produced at the early stage of 
relativistic heavy ion collisions. At the high collision en- 
ergies RHIC and top SPS, the yields of different hadron 
species have shown an explicit "quark-level" correlation 
in quark (re-)combination/coalescence scenario ifTsl - flTll . 
In this paper, we make an energy scan from RHIC ener- 
gies vtwv = 200, 130, 62.4 GeV to SPS energies = 
158,80,40,30,20 AGeV to study at which collision en- 
ergy this "quark-level" correlation of hadron yields first 
breaks. In particular, we define two correlation quantities 

A = xipr= and B - 4^^-§-, sensitive to quark degrees 
of freedom. The values of A and B are equal to one in 
the framework of quark (re-)combination/coalesce, inde- 
pendent of models. The deviation of A and B from one or 
not can be regarded as a possible signal of deconfinement 
in heavy ion collisions. We apply a quark combination 
model, which can exclusively describe hadron production 
and well reproduce the yields and momentum spectra of 
final-state hadrons in relativistic heavy ion collisions ifisl- 
I23I1 . to carry out the concrete calculations. 

The paper is organized as follows. In Sec. II, we present 
the relations among yields of various hadrons in quark (re- 
)comibination/coalescence scenario. In Sec. Ill, we use the 
quark combination model to calculate the yields of various 
hadrons, their yield ratios and correlation quantities A = 

^jT-= and B - ^ ^_ at midrapidity in the most central 
A+A collisions at different energies. Sec.IV summaries 
our work. 



II. HADRON YIELDS IN QUARK 
(RE-)COMBINATION/COALESCENCE SCENARIO 



'Electronic address: [shaofl@mail.sdu.edu.cnl 



Let us start from the general inclusive formula of 
hadron production in quark (re-)combination/coalescence 
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Nm 
Nb 



= J dpidp2Fg^ipi,p2)'RM(PuP2) (1) 
= J dpidp2dp3Fggg(puP2,P3)'RB(Pl,P2,P3)i2) 



Here Fgg (Fqgq) is the joint quark-antiquark (three quark) 
distribution. CRg) is the combination function which 
stands for the formation probabiUty of quark antiquark 
(three quarks) into a meson (baryon), dominated by chro- 
modynamics. In sudden approximation, it is equal to the 
overlap between two (three) quark wave functions and 
the wave function of meson (baryon). Neglecting exotic 
(multi-quark) states, mesons and baryons exhaust all fate 
of quarks and antiquarks. One reaches the following rela- 
tions : 2 A?M + 3 2 A^B = 2 Nq and 2 A^m + 3 2 A^^ = 2 Nq, 
where A^^ (A^^^) is the quark (antiquark) number of flavor 
qiq). Extracting A^^^ and A^^ from the joint two (three) quark 
distribution F^^ (Fqqq) and putting them out of the integral, 
one has the following schematic relations between hadron 
yields and quark numbers after integrating over quark mo- 
menta 



NM(qq) « CMNqNq, 



NBUm) <^ CeNqNqNq, (3) 



and quark number conservation will fix the proportionality 
coefficient. The effects of combination function on hadron 
yields are characterized with the factors Cm and Cb, and 
Cm - Cjj and Cb - C-g are assumed. 

One realized method of quark number conservation dur- 
ing combination is adding a factor bq for each quark flavor 
in above equations, as did in ALCOR model ll24ll . 

A^M(w) = CM(bqNq){bqNq), (4) 
NB(,qq) = CB(bqNq){bqNq)(bqNq). (5) 

Then bq can be uniquely determined by quark number con- 
servation. According to Eqs.lHli and (|5]), we have the fol- 
lowing relations between hadrons and the corresponding 
antihadrons 



P ^ bqq^^ 
P bqq 

A ^ ^bqq^^bjs 
A bqq bss' 



bqq bsS 

_ ^bss^2 bgl 
bjs bqq' 



Here, we use particle symbols stand for their numbers for 
short, q for light quark number and s for strange quark 
number Hiding the quark content in hadron yield, we 
obtain the following interesting relations among different 
hadron species 
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We define correlation quantities A and B as follows 



A = 



Ar 



Ar p 



B 
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If the quark matter exists and hadronizes via sudden (re- 
)combination/coalesce in heavy ion collisions, A and B 
should be equal to one for directly produced hadrons. It 
is a general result under the constraint of quark number 
conservation, which is independent of specific models. 

Different from the well-known recombination and co- 
alescence models ifisl [Till , the quark combination model 
(H, Iz^ is unique for its combination rule. The main idea 
of the combination rule is to line up the (anti)quarks in 
a one-dimensional order in phase space, e.g., in rapidity, 
and then let them combine into initial hadrons one by one 
according to this order ifisl . Three (anti)quarks or a quark- 
antiquark pair in the neighborhood form a (anti)baryon or 
a meson, respectively. At last all quarks and antiquarks 
are combined into hadrons. The relations between hadron 
yields and the corresponding quark numbers are easily ob- 
tained lfl9ll . With this rule, the model can give the yields 
and momentum distributions of all hadrons (included in the 
model) in an event, possessing some exclusive nature. The 
decay of the short-life resonances is systematically taken 
into account in the model. The model has been realized in 
Monte Carlo program and has described many properties 
of hadron production in relativistic heavy ion collisions 



III. RESULTS AND DISCUSSIONS 

In this section, we use the quark combination model to 
calculate the hadron yields and their ratios as well as cor- 

^'"^ andB = 



at midra- 



(6) 



relation quantities A = ^, Ak+ 

pidity in the most central A+A collisions at yfsjm - 
200, 130, 62.4 GeV and Ebeam = 158, 80, 40, 30, 20 
AGeV. The predictions at LHC are also presented. The 
necessary input of the model, i.e., quark distribution just 
before hadronization, can be obtained by applying the hy- 
drodynamics to describe the evolution of hot and dense 
quark matter just before hadronization, see Appendix |A] 
for details. 



Hadron yields, their ratios and correlation quantities A 
andB 



Table U shows the hadron density dN/dy at midrapidity 
in the most central A+A collisions at RHIC energies 200, 
130, 62.4 GeV and SPS energies 158, 80, 4(1 30, 20 AGeV. 
The experimental data are taken from Refsll9l l28l - (38ll . The 
agreement between the calculation results and the experi- 
mental data is good except at 20 AGeV where y^/nd f is 
far greater than one. We note that ALCOR model ll24ll 
can also describe well the hadron yields at top SPS energy, 
which gives a cross-verification of quark combination. On 
the contrary, HIJING or HIJING/B model which is in com- 
pliance with the hadronic scenario for the early evolution 
of heavy ion collisions via a fragmentation hadronization, 
can not self-consistently explain the data of multi-strange 
hadrons ||39i|40|l. 
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TABLE I: The calculated hadron yields dN/dy at midrapidity in the most central A+A collisions at different energies. The experimental 
data are taken from Refs(^[2^^. 





Au+Au 200 GeV 


Au+Au 130 GeV 


Au+Au 62.4 GeV 


Pb+Pb 158 AGeV 




data 


model 


data 


model 


data 


model 


data 


model 


+ 

n 


286.4 ± 24.2 


281.0 


276 ± 3 ± 35.9 


267.2 


233 ± 17 


227.4 


170.1 ±0.7 ±9 


165.1 




281.8 + 22.8 


281.8 


270 ±3.5 ±35.1 


270.5 


237 ± 17 


233.5 


175.4 ±0.7 ±9 


176.5 




48.9 ± 6.3 


48.6 


46.7 ± 1.5 ±7.0 


45.2 


37.6 ±2.7 


38.3 


29.6 ± 0.3 ± 1.5 


27.2 




45.7 ± 5.2 


46.1 


40.5 ±2.3 ±6.1 


42.4 


32.4 ± 2.3 


32.2 


16.8 ±0.2 ±0.8 


16.2 


P 


18.4 ±2.6 


17.0 


28.7 ±0.9 ±4.0 


25.7 


29.0 ±3.8 


29.1 


29.6 ± 0.9 ± 2.9 


30.1 


P 


13.5 ± 1.8 


12.5 


20.1 ± 1.0 ±2.8 


18.2 


13.6 ± 1.7 


13.5 


1.66 ±0.17 ±0.16 


1.91 


A 


16.7 ±0.2 ± 1.1 


15.3 


17.3 ± 1.8 ±2.8 


14.5 


14.9 ± 0.2 ± 1.49 


13.7 


10.9 ± 1.0 ± 1.3 


13.3 


A 


12.7 ±0.2 ±0.9 


12.1 


12.7 ±1.8 ±2.0 


10.9 


8.02 ±0.11 ±0.8 


7.28 


1.62 ±0.16 ±0.2 


1.66 




2.17 ±0.06 ±0.19 


2.05 


2.04 ±0.14 ±0.2 


1.93 


1.64 ±0.03 ±0.014 


1.67 


1.44 ±0.10 ±0.15 


1.18 




1.83 ± 0.05 ± 0.20 


1.69 


1.74 ±0.12 ±0.17 


1.53 


0.989 ± 0.057 ± 0.057 


1.02 


0.31 ±0.03 ±0.03 


0.23 


n- 


0.53 ± 0.04 ± 0.04 


0.56 


0.56 ±0.11 ±0.06 


0.52 


0.356 ±0.046 ±0.014 


0.379 


0.14 ±0.03 ±0.01 
0.07 ± 0.02 ± 0.01 


0.11 
0.04 


X^lndf 


2.8/8 


2.0/8 


2.1/8 


8.5/9 





Pb+Pb 80 AGeV 


Pb+Pb 40 AGeV 


Pb+Pb 30 AGeV 


Pb+Pb 20 AGeV 




data 


model 


data 


model 


data 


model 


data 


model 




132.0 ±0.5 ±7 


129.9 


96.6 ± 0.4 ± 6 


97.8 


83.0 ± 0.4 ± 4.2 


84.8 


72.9 ±0.3 ±3.6 


73.6 


n" 


140.4 ± 0.5 ± 7 


141.8 


106.1 ±0.4 ±6 


110.1 


96.5 ±0.5 ±4.8 


99.3 


84.8 ± 0.4 ± 4.2 


85.0 


k* 


24.6 ± 0.2 ± 1.2 


23.8 


20.1 ±0.3 ± 1.0 


19.4 


21.2±0.8!i:^ 


21.6 


16.4 ±0.6 ±0.4 


16.6 


k- 


11.7 ±0.1 ±0.6 


11.9 


7.58 ±0.12 ±0.4 


7.27 


7.8 ±0.1 ±0.2 


7.3 


5.58 ±0.07 ±0.11 


5.01 


P 


30.1 ± 1.0 ±3.0 


31.5 


41.3 ± 1.1 ±4.1 


38.6 


42.1 ±2.0 ±4.2 


37.7 


46.1 ± 2.1 ±4.6 


37.0 


P 


0.87 ± 0.07 ± 0.09 


0.83 


0.32 ± 0.03 ± 0.03 


0.36 


0.16 ±0.02 ±0.02 


0.17 


0.06 ± 0.01 ± 0.006 


0.07 


A 


13.5 ±0.7 ±1.0 


14.5 


15.3 ± 0.6 ± 1.0 


14.7 


14.7 ± 0.2 ± 1.2 


14.2 


13.4±0.1 ± 1.1 


11.5 


A 


1.06 ±0.08 ±0.1 


0.95 


0.42 ± 0.04 ± 0.04 


0.49 


0.21 ±0.02 ±0.02 


0.20 


0.10 ±0.02 ±0.01 


0.08 




1.22 ±0.14 ±0.13 


1.23 


1.15±0.11±0.13 


1.10 


1.17±0.13±0.13 


1.38 


0.93 ±0.13 ±0.10 


0.92 




0.21 ± 0.03 ± 0.02 


0.15 


0.07 ±0.01 ±0.01 


0.08 


0.05 ± 0.01 ± 0.01 


0.07 




0.02 


fl" + f2 




0.13 


0.10 ±0.02 ±0.02 


0.09 




0.14 




0.07 


x'/ndf 


2.9/7 


3.0/8 


5.6/7 


15.3/6 



Fig.[T] shows the ratios of antihadrons to hadrons at 
midrapidity as the function of collision energy. The filled 
triangles are the calculated results, and the experimental 
data are taken from Refs[9, 28,-30, 32-34, 36,-38i|. These 
ratios are mainly influenced by the no-zero baryon num- 
ber density. As the collision energy increases, the nuclear 
transparency power becomes strong and baryon number 
density at midrapidity becomes small. This results in a 
rapid increase for ratios of K^/K^, p/ p, A/A with the in- 
creasing collision energy, tt^/tt^ follows a different pat- 
tern. At low collision energies n^/n^ is slightly higher 
than one, while it is close to one at high energies. This is 
caused by the asymmetry of decay contribution from hy- 
perons and anti-hyperons (e.g. A p tt ). As the colli- 
sion energy increases, the yields of hyperons are close to 
that of anti-hyperons and their decay contributions to pion 
yields are almost the same and therefore the ratio of /n^ 
is close to one. 

Fig. |2] shows the ratios of baryons to mesons p/n and 
the strangeness ratios K/tt at midrapidity as the function 
of collision energy. The filled triangles (up and down) are 



the computed results, and the experimental data are taken 
from RefsH [M [US The big splits between 

p/n^ and p/n^ and between jii^ and K^/tt^ at low col- 
lision energies is due to the high baryon number density. 
At 20 AGeV, the result of p/n^ deviates seriously from 
the data. This is probably because the participant nucleons 
are not broken completely in collisions. These nucleon 
fragments deposited in midrapidity region lead to the ex- 
tra contribution of proton production besides those from 
quark combination and lead to the excessively high ratio 
of p/tt^. 

Following the experimental data in Table H] A and B at 
different collision energies are evaluated and the results are 
presented with squares in Fig|3] We find that the data of 
A and B at RHIC energies are almost equal to one while 
at SPS energies they deviate from one; particularly, at 20 
AGeV the value of A amounts to two, seriously deviating 
from one (the data for B are unavailable). This is proba- 
bly because the decay of resonances will blur A and B to 
a certain degree. In order to explore the decay effect, we 
use the quark combination model to compute the values 
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FIG. 1 : (Color online)The yield ratios of antiliadrons to liadrons 
at midrapidity as the function of the collision energy. The filled 
symbols are the calculated results, and the experimental data, 
open symbols with error bar, are from Refsl3. l28l43ft[32l434l . [3a - 
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FIG. 2: (Color online)The values of the relative production yields 
p/n and k/n at different energies. The triangles down are the 
numerical results for p/n* and k* In* , and the triangles up are 
the numerical results foip/n^ and k^ jn^ . The experimental data, 
open symbols with error bar, are from Refs[9. 28 . 32 - 34 . 36„.37ll . 



of A and B for the directly produced hadrons and the final- 
state hadrons, respectively. The dashed lines are the results 
of directly produced hadrons. Just as analyzed above, the 
dashed lines keep an invariant value of one for both A and 
B, independent of collision energy. The very small fluctua- 
tions of dashed lines are due to the fact that the rapidity dis- 
tributions of the formed hadrons are slightly different from 
those of quarks, which leads to a small amount of hadrons 
formed by midrapidity quarks escape from the midrapid- 
ity region ifisll . The filled triangles down are model re- 
sults of final-state hadrons. We see that as collision energy 
is greater than 30 AGeV these triangles agree well with 
the experimental data within statistical uncertainties. Re- 
moving the part of resonance decay, the data of A and B 
for directly produced hadrons should be equal to one, re- 
spectively. This suggests the existence of the "quark level" 
correlation of directly produced hadrons at these collision 



energies. However, at 20 AGeV the value of A calculated 
via quark combination for final state hadrons seriously de- 
viates the data. This in itself indicates the quark degrees of 
freedom do not represent a decisive factor in hadron pro- 
duction at 20 AGeV. 
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FIG. 3: (Color online)The correlation quantities A and B as the 
function of collision energy. The experimental data, open sym- 
bols with error bar, are from Ref s L 9^. . 2 8.- 3 0. . 3 2-, 3 8i] . 

From the above analysis of hadron yields, hadron ratios 
and correlation quantities A and B, we find that as the col- 
lision energy is greater than or equal to 30 AGeV the quark 
combination can describe reasonably well all these observ- 
ables which indicates the existence of constituent quark de- 
grees of freedom in the energy region. When the collision 
energy drops to 20 AGeV, however, the quark combination 
mechanism can not self-consistently describe these quan- 
tities, in particular the correlation quantity A. Note that in 
order to best fit the data of hadron yields, we have adjusted 
the strangeness parameter As (see Table llVt to the quite 
high values at 20, 30 AGeV which are far greater than 
the saturated values at higher collision energies and are 
also greater than those at lower energies ll4lll . But even so 
the quark combination can not self-consistently explain the 
data at 20 AGeV. In addition, the high strangeness entan- 
gled with baryon density leads to the Sawtooth-like shape 
for calculated final-state A and B at low SPS energies. In 
fact this peak behavior of strangeness (i.e. K* In* ratio or 
As determined mainly by the former ) at low SPS energies 
has been interpreted in Ref 01 as a result of onset of 
deconfinement, i.e. the result of strangeness carrier chang- 
ing from strange hadrons to strange quarks at the onset of 
deconfinement. The failure of quark combination energy 
indicates the (partial) disappearance of constituent quark 
degrees of freedom and closely relates to the onset of de- 
confinement around 30 AGeV observed by NA49 Collab- 
oration i^. 



B. Predictions at LHC 

With the increasing collision energy, the strangeness As 
of the hot and dense quark matter in Table |IV] tends to 
be saturated and the squared sound velocity approaches 
an ideal value 1 /3 and the baryon number density «o de- 
creases regularly. So we take As-QA1>, — 1/3 and 
«o = O.OfmT^ at yfsNN - 5.5 TeV The initial entropy den- 
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sity of the central point at the beginning of hydrodynamic 
evolution is taken to be 256.5 /fm^ according to Eq. ( |A7t . 
In Tables HI] and Hill we present our predictions of hadron 
yields and hadron ratios at LHC energy ( -y^siviv = 5.5 TeV). 
Our computed particle ratios are consistent with the results 
of the thermal model in RefH. 
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Appendix A: Quark distribution just before hadronization 



TABLE II: The predicted dN/dy of identified hadrons at midra- 
pidity in the most central (0 - 5%) Pb+Pb collisions at = 
5.5 TeV. 



n 



E- Q- 



Q. 



435.5 435.5 71.6 71.6 34.6 34.6 20.5 20.5 2.9 2.9 0.4 0.4 



TABLE III: The predicted hadron ratios at midrapidity in the 
most central Pb+Pb collisions at s/s^n = 5.5 TeV. 

p/n+ pin- k+/7i+ k-/n- A/n- E-fn- 0.-/71- 
0.079 0.079 0.165 0.165 0.047 0.007 0.001 



IV. SUMMARY 

In this paper, we used the quark combination model to 
make a systematical study of the hadron yields in heavy 
ion collision in a broad collision energy region. At the col- 
lision energies where the deconfined quark matter has been 
created, the yields of various hadrons after hadronization 
have some correlations inherited from the early quark de- 
grees of freedom. We investigate at which collision energy 
this "quark-level" correlation of hadron yields first breaks. 
We apply the hydrodynamics to describe the evolution of 
deconfined quark matter and to obtain the quark distribu- 
tion just before hadronization, then utilize the quark com- 
bination model to describe hadronization. We find that as 
the collision energy is greater than or equal to 30 AGeV, 
the quark combination well reproduce the yields of vari- 
ous hadrons, their ratios and correlation quantities A and 
B; however, as the collision energy drops to 20 AGeV, the 
mechanism can not self-consistently describe these quan- 
tities. This indicates that the constituent quark degrees 
of freedom represent a decisive factor in thermal hadron 
production above 30 AGeV and seem to be invalid at 20 
AGeV. It is related to the onset of deconfinement observed 
at collision energy 20-30 AGeV. Finally, we predict the 
yields of various hadrons and their ratios in the most cen- 
tral Pb+Pb collisions at -Jsjm = 5.5 TeV. 
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The quark distribution just before hadronization is 
needed for the quark combination to describe the hadron 
production in relativistic heavy ion collisions. Here, we 
use relativistic hydrodynamics ll43ij45ll to describe the 
time-space evolution of the hot and dense quark matter 
before hadronization. The evolution equation of hydro- 
dynamics follows from the local conservation laws for en- 
ergy, momentum, and other conserved charges, e.g. baryon 
number. 



5„r''''(x) = 0, (v = 0,1,2,3) 



(Al) 



(A2) 



d,nx) = 0, 
by inserting the ideal fluid decomposition 

Tf^ix) = ie{x) + p{x)y^'{x)u\x) - g^^pix), (A3) 



fix) = n{x)u''(x). 



(A4) 



Here, u^^ix) - yil. 



,V;) withy = 



is the local four velocity of a thermalized fluid cell; e{x) 
is the energy density, p(x) the pressure, and n{x) the con- 
served number density. 

As the energy density of the fluid cell drops to 1.0 
GeV/fm^* (a common criteria of phase transition from Lat- 
tice OCD ll46ll ). we stop the hydrodynamic evolution and 
let the constituent quarks and aritiquarks freeze out accord- 
ing to Cooper-Frye formalism ||47I1 



^ ( flip ■ uix), x)p ■ d^o-ix). 



dNi dNi 

d^p dyprdprdip (2: 

(A5) 

where d^cr{x) is the outward normal vector on the freeze- 
out surface 2(jc), gi the degeneracy factor of quarks {gi = 
6). The phase-space distribution / in the formula is taken 
to be a local equilibrium distribution. 



fi{E,x) 



1 



exp[(£-//,W)/r(x)] + l 



(A6) 



We consider only the freeze-out of light and strange quarks 
and antiquarks. The three chemical potentials ^q, fig and 
fis - /^j at freeze-out can be determined uniquely by the 
global conservation of energy and baryon number plus an 
ancillary constraint of strangeness. 

We study in this paper the hadron production in midra- 
pidity region only. The code of Kolb ll48l - l50ll for 2+1- di- 
mensional hydrodynamics with the longitudinal boost in- 
variance is used to simulate the evolution of quark system 
before hadronization. 
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TABLE IV: The values of the initial baryon density ng, strangeness factor /I, and squared sound velocity in central A+A collisions 
at different energies. 



Energy 200GeV 


130GeV 


62.4GeV 


158AGeV 


SOAGeV 40AGeV 


30AGeV 20AGeV 


no (fm-^) 0.30 


0.34 


0.74 


1.31 


1.46 1.64 


1.74 1.56 


A, 0.43 


0.43 


0.43 


0.44 


0.50 0.57 


0.80 0.70 


c? 1/3.1 


1/3.4 


1/4.0 


1/6.0 


1/6.0 1/6.0 


1/6.0 1/6.0 



These frozen-out quarks and antiquarks with momen- 
tum distributions in Eq. iA5i are hadronized by the quark 
combination model 11 15l 12511 . In this method of hydro- 
dynamics -I- quark combination, the quark combination 
model is default setting while some inputs for hydrody- 
namics need fixing. The first is initial entropy density sq 
of the central point at the beginning of hydrodynamic evo- 
lution. The distribution of transverse entropy density is 
determined by optical Glauber modell^sSl. We apply the 
following empirical formula of i()|27!] as the function of 
collision energy yfs^N 

soi = 6.99 X 10^^^(312.5 logi,, - 64.8)^^1 

(A7) 

The second is the initial baryon number density no at cen- 



tral point fixed by the data of the net-proton rapidity den- 
sity. The distribution of transverse baryon number den- 
sity is also determined by optical Glauber model. The 
third is the equation of state p - p{e) for the hot and 
dense quark matter which is taken to be the simplest pat- 
tern p - c^e. The squared sound velocity c] is obtained 
by fitting the data of transverse momentum spectra of pro- 
tons. The fourth is the strangeness of hot and dense quark 
matter denoted by the factor /Ij = 2{ss)l{uu + dd). In order 
to make a better description of strange hadrons, we regard 
it as a free parameter in the present paper. The values of 
«o, and Cj in central A+A collisions at different ener- 
gies are shown in Table II VI We note that the extracted As 
around 30 AGeV are quite high which are consistent with 
the analytic results of thermal model iiHl . 
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